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The connection between metallicity and metal-line kinematics in 
(sub-)damped Lyman-or systems 
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ABSTRACT 

A correlation between the metallicity, [M/H], and rest-frame Mgii equivalent width, W? 2196 , 
is found from 49 DLAs and strong sub-DLAs drawn from the literature over the redshift 
range 0.2 < z a t, s < 2.6. The correlation is significant at 4.2 cr and improves to 4.7 cr when 
the mild evolution of [M/H] with redshift is taken into account. Even when including only 
the 26 DLAs (i.e. exlcuding sub-DLAs) which have Zn-metallicities and W? 2196 > 0.7 A, the 
correlation remains at > 3 cr significance. Since the Mg n /12796 transition is predominantly 
saturated in DLAs (which always have W} 2196 greater than 0.3 A), W} 2196 is far more sensitive 
to the kinematic spread of the metal velocity components across the absorption profile than 
it is to [M/H]. Thus, the observed [M/HJ-W^ 2796 correlation points to a strong link between 
the absorber metallicity and the mechanism for producing and dispersing the velocity com- 
ponents. We also note that approximately half of the 13 known H2 absorbers have very high 
Wj )27% and very broad velocity structures which show characteristics usually associated with 
outflows. Follow-up ultraviolet- and blue-sensitive high-resolution spectra of high- W^ 2796 sys- 
tems, initially identified in low-resolution spectra, may therefore yield a large number of new 
H2 discoveries. 
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1 INTRODUCTION 

The damped Lyman-cf absorbers (DLAs) observed along sight- 
lines to quasars (QSOs) are potentially important pieces in the 
puzzle of galaxy formation. They are the highest column density 
absorption systems, with iV(Hi) > 2 x 10 20 cnr 2 , and the first 
DLA surveys quickly established that they contain a large frac- 
tion of the high-z baryons a vailable for star formation { Wolfe et al.l 
ll98d : lLanzetta et alJI 199lh . The early estimates of the neutral gas 
mass density at absorption redshifts z a bs ~ 3 were similar to the 
estimates of the mass density in stars at z = 0. It was there- 
fore natural to assume that DLAs act as the neutral gas reser- 
voirs for star formation. However, firmly establishing the link be- 
tween DLAs and star formation relies on our understanding of the 
absorber-galaxy connection. Direct imaging of DLAs (or strong 
Mgii systems; see below) at z a bs < L5 reveals the hosts to be a 
mix of i rregulars, dwarfs, spirals and low surface-brightness galax- 
ies (e g. iBergeron & Bpisselll99ll : iLe Brun et al]|l997l : iRao et all 
120031 : IChen & La nzetta 2003]). Th is is further borne out by a blind 
21 -cm emission survey at z = ( Rya n- Weber et al. 1 120031) . How- 
ever, quantifying the mix of DLA host-galaxy morphologies from 
low- to high-redshift is fraught with selection effects, not least 
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of which are luminosity bias against faint galaxies and the prox- 
imity of some host-galaxies to the QSO sight-line. Indeed, very 
few DLA host-galaxies have been identified at z a bs <: 2 (e.g. 



M0ller & Warre: 
2005tlWeatherlevet 



1993 1; iDiorgovski et alj fl996l : IChristensen et all 
alj|2005h 



The Mgii /12796/2803 doublet has proved very important 
for exploring the absorber-galaxy connection at 0.3 < z a bs 
2.6 where it is observable with optical QSO absorption spec- 
troscopy. Indeed, many of the z abs < 1.5 imaging studies above ini- 
tially identified absorpti on systems via the easily recognised Mg 11 



doublet (see also, e.g., Lanzetta & Bowenl 19901: [Bergeron etal 



1992l:ISteidel & Sargendl992l:lDrinkwater et alJl993; Steid el et al 



1994h . Also, since the ionization potential of Mg 1 is < 13.6 eV but 
that of Mg 11 is > 13.6 eV, the latter traces cold gas. Attempts have 
therefore been made to id entify DLA candidates at low-z abs vi a 
their Mgn absorption (e.g. IRao et al]|l995L 120061 ; lEUisonl l2006). 
The Mgn doublet transitions are also among the strongest metal 
lines observed in QSO absorbers, meaning that (i) the Mg n /12796 
rest-frame equivalent width in DLAs is typically W; t27% > 0.3 A, 
allowing such systems to be identified in relatively low-resolution 
spectra, and (ii) most velocity components across the Mg 11 profile 
are saturated in DLAs. This implies that W r ' 127% is mainly sensitive 
to the kinematic extent of the velocity components, Av, rather than 
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their column density. A good demonstration of this is figure 3 of 
lEllisodJ2006h . 

In this paper we search for a relationship between W? 2196 
and the metallicitjQ, [M/H], of DLAs and strong sub-DLAs. Since 
W; 2796 is a measure of the kinematic spread of the cold, metal- 
bearing velocity components, such a relationship may constrain 
the mechanisms for producing and dispersing the metals. This may 
guide our understanding of the physical nature of DLAs. 

There currently exist only indirect and possibly contradic- 
tory constr aints on a [M/HI -IV; 12796 relationship for DLAs and 
sub-DLAs. iNestor et ail d2003h suggest that [M/H] increases with 
increasing W/ 2796 . They identified strong Mgn systems in Sloan 
Digital Sky Survey (SDSS) QSO spectra over the redshift range 
0.9 < Zabs < 2.0 and constructed composite spectra in the ab- 
sorber rest-frame. The Znn lines were stronger in the composite 
of W/ 2796 > 1.3 A absorbers compared with those in the composite 
of 1.0 < W; 12796 < 1.3 A systems. Assuming that the mean N(Hi) 
is not significantly different in these two regimes, they conclude 
that some e vidence exists for a [Zn/H]-lV; t2796 correlation. On the 
other hand, lYork et al.1 d2006h find less direct evidence for an anti- 
correlation between [Zn/H] and VV; 12796 . From a series of composite 
SDSS spectra with various mean W; 12796 they find that the dust- 
reddening, E(B - V), caused by the absorbers rises sharply with 
increasing W* 2196 . While the Znn line-strengths also increase, they 
show a slower increase than E(B- V) over the same range of W r ' 12796 . 
The composite spectra also suggest a dust extinc tion curve t ypical 
of the Small Magellanic Cloud (SMC), leading lYork et all to as- 
sume a constant dust-to-gas ratio, i.e. A'(Hi) oc E(B - V), as ob- 
served in the SMC. Thus, since the Znn column density increases 
slow er than the iV(H i) inferred from E(B - V) under this assump- 
tion, lYork etalj 12006) argue that [Zn/H] decreases with increasing 
W r ' 12796 . We discuss these results further in light of our new results 
in Section[3] 

This paper is organised as follows. Section[2]explains our sam- 
ple selection and explores an evident [M/HJ-W,; 12796 relationship 
in that sample. Section [3] compares our results w i th the [M/H]- 
Av correlation recently discovered bv lLedoux et al.l d2006h at (gen- 
erally) higher redshifts. We also discuss whether our observed 
[M/HJ-W^ 2796 relationship can be taken as evidence for a corre- 
lation or anti-correlation between absorber metallicity and host- 
galaxy mass. Finally, we discuss whether a large fraction of strong 
Mgn absorbers might arise in outflows from relatively low mass 
galaxies. In particular, about half of those absorbers in which H 2 
has been detected show some evidence of an outflow origin. 



2 CORRELATION BETWEEN METALLICITY AND 

MG II EQUIVALENT WIDTH IN DLAS AND SUB-DLAS 

2.1 Sample definition 

We identified DLAs [N(R i) > 2 x 10 20 cnr 2 ] and strong sub-DLAs 
[A'(Hi) > 3 x 10 19 c m -2 ] with existing me tallicity measurements 
in the compilation of Kulk arni et alj |2005) and several other, typ- 
ically more recent, sources (as cited in Table [TJ. For these systems 
we searched the literature for W? 2796 measurements; Table Q] lists 
those systems for which we found both [M/H] and W? 2196 . Un- 
certainties for some W? 2796 values could not be found in the lit- 



1 Metallicity is defined as the heavy element abundance with respect 
to hydrogen, relative to that of the solar neighbourhood: [M/H] = 
log[iV(M)/W(H)] - log[iV(M)/iV(H)] 



erature (see Section [2~2\ . For systems without a published W; 12796 
value we searched for publicly available SDSS (R ~ 2000) or 
Very Large Telescope Ultraviolet and Visual Echelle Spectrograph 
(VLT/UVES; R > 40000) spectra from which an IV/ 2796 measure- 
ment is easily made. For SDSS s pectra, our method for obtaining 
Wf 196 is described in detail in lBouche et"ai1 ( 120061) . UVES spec- 
tra were obtained from the ESO data archive and reduced using a 
version of the UVES pipeline which we modified to improve the 
flux extraction and wavelength calibration. The extracted echelle 
orders from all exposures were combined, using uves popleiQ, with 
inverse-variance weighting and a cosmic ray rejection algorithm, 
to form a single spectrum with a dispersion of 2.5kms~' pixel -1 . 
^12796 was men measurec j by direct summation of pixels across the 
Mgn /12796 profile. The uncertainties in W? 2196 quoted in Table Q] 
for these UVES spectra are generally dominated by uncertainties 
in the continuum level around the Mg n /12796 line. In some cases 
where Wf 2196 is so high that the /12796 profile is blended with that 
of the /12803 line (typically VV; 12796 > 3.5 A), larger uncertainties 
are reported. 

For most of the 49 absorbers in Table [T] Zn-metallicities were 
available in the literature. In 13 cases we have used the Fe or Cr 
abundance in the absence of a Zn measurement. It is well known 
that t hese elements are depleted onto dust grains (e.g. lPettini et al.l 
1 19971) and so we corrected the Fe and Cr abundances with the sim- 
ple prescription of IProchaska et alj d2003al) : [M/H] = [Fe/H] + 
0.4 = [Cr/H] + 0.2. In Section[2~2]we test the possible impact this 
might have on our results. Most of the metallic i ties in Table[T]have 
been corrected to the solar scale of Lodders (2003). The values 
quoted in the literature were used for Q0235+164, Q 1622+238, 
Q 0458-020, Q 1122-1649, Q 1157+014 since we could not iden- 
tify which solar values the original authors used. Since the correc- 
tions are typically < 0. 1 dex, these 5 exceptions have negligible 
impact on our main results. 



2.2 The Wf 796 -[M/H] correlation 

Figure[T]left) shows all the [M/HJ-VV; 12796 pairs. TableQ]gives the 
formal statistical errors on both quantities as reported in the litera- 
ture. In cases where no error-bars could be found in the literature, 
uncertainties of 0. 10 dex and 0.05 A were assumed as these typified 
the rest of the sample. For Fig.QJleft) we added 0.05 dex and 0.03 A 
in quadrature to the initial [M/H] and W? 2196 uncertainties to take 
into account likely continuum fitting errors which are not included 
by some authors or, indeed, by our own H^ 2796 measurements in 
SDSS spectra. We regard the points in Table[T]with these increased 
errors as our fiducial sample. Kendall's r for the [M/HJ-VV^ 2796 cor- 
relation evident inFig.Qjhas aprobability P(t) = 2.6xl0~ 5 of being 
due to chance alone; assuming Gaussian statistics, the correlation 
is significant at >4.2-<r. The correlation is robust to the removal of 
distinct sub-samples, some of which are discussed below, as shown 
in Table|2] 

The absorbers cover quite a large redshift range, 0.2 < z a bs < 
2.6, so evolution in DLA metallicities may have somewhat washed 
out the correlation we observe. For example, if absorbers at lower 
z a bs are, on average, more metal-rich than at higher z a b S! then this 
would cause additional dispersion in [M/H] at all values of W; ,27% . 
This would at least decrease the statistical significance of the ob- 
served correlation and, depending on the redshift distribution of the 
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Metallicity and kinematics of (sub-)DLAs 



Table 1. DLAs and sub-DLAs where both Mgn ,12796 A equivalent width, W? 2196 , and metallicity, [M/H], have been measured. B1950 QSO 
names are given. A'(Hi) is the total neutral hydrogen column density [cm -2 ], z aDS is the DLA redshift and z cm is the redshift of the background 
QSO. Errors on W^ 2796 and [M/H] are quoted from the liter ature but in many cases do not include continuum-fitting uncertainties. [M/H] 
is measured relative to the solar abundances of Lodders 1 2003) in most cases (see text). The element used to d etermine [M/H] is Zn w here 
possible; when this is not available we corrected for dust depletion via [M/H] = [Fe/H] + 0.4 = [Cr/H] + 0.2 (see lProchaska et alj2003ah . 

QSO logW(Hi) Zabs z cm W, 2796 [A] Ref. [M/H] M Ref. Notes 
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2.066 


3.015 


1.80 


E06b 


-0.86 ±0.10 


Zn 


D04 


e 


2343+125 


20.3 


2.4310 


2.549 


3.22 ±0.15 


UVES 6,12,17 


-0.74 ± 0.08 


Zn 


D04 


a, e 


2348-1444 


20.6 


2.279 


2.94 


0.32 


E06b 


-1.84 


Fe 


D06 


d 



Notes: a H2 absorption detected; ''21-cm absorption detected; c 21-cm absorption not detected; ''radio-loud but not searched for 21-cm absorption; e radio-quiet. '60.A-9022 UVES 
commissioning; 2 65.O-0063 Ledoux; 3 65.0-O158 Pettini; 4 65 P-0038 Srianand; 5 66.A-0212 Reimers; 6 66.A-0594 Molaro; 7 67.A-0O22 D'Odorico; S 67.A-0146 Vladilo; 9 68 A-0106 
Petitjean; 10 68.A-0361 Lopez; n 68.A-0600 Ledoux; 12 69.A-0204 D'Odorico; 13 69.B-0108 Srianand; 14 7O.A-0O17 Petitjean; 15 70.B-0258 Dessauges-Zavadsky; 16 072.A-0100 Murphy; 
17 072.A-0346 Ledoux; 18 073.B-0787 Dessauges-Zavadsky; 19 074.A-0201 Srianand. 

References: SM79 ISpinrad & McKeJl979tl . B87 iBoulade et alll987h . W93 IWolfe et all 1981 . LB92 lLanzetta & Bowenl 1992fc . MY92 iMeyer & Yorkll992t . SS92 ISteidel & Sargenj 
119921 ). A94 lAldcroftetall!994 . M95 iMever et alll995t . P97 iPettM et alll997t . B98 iBpisse et alll998l) . L99 lLopez et all 19991 ), P99 iPettiniet all 19991 ). PW99 iProchaska & Wolfd 
1 19991) . COO I Chur chill et al. 2000), D00 i de la V area et allEoQcl) . POO IPettini et allboOfj) . RT00 iRao & Tuntshekll200fl) . P01 IProchaska etaj]l200 lb . P02 IPetitiean et ai]|2002l) . L02a 
Ledoux et al. 2002). L02b Ledoux et al. 2002), C03a Centurion et al. 2003), C03b i Churchill et al. 2003), D03 . Dessauges-Zavadsky et al. 2003). L03 Ledoux et al. 2003). LE03 
lLopez & Ellisorll2003l) . P03a [Sochaska et al. 2003al), P03b [Prochaska et al. 2003tl). R03 iRyabinkov et al. 20(3) [and references therein], D04 I Dess auges-Zavadsky et al. 20041) , J04 
lJunkkarin en et al. 2004). K04 i Khare et al. 2004), P04 ■ Peroux el all2004 . A05 [Xkerman et al. 2 0o3) . K05 [Kujkaini el al. 2005). L05 i Lopez et al. 2005), D06 iDessauges-Zavadskv et all 
120061) . E06a lEllisorl2006l) . E06b (S. E. Ellison, private communication), P06a IPeroux et al. 20o3), P06b IProc haska et al. 2006). R06 i Rao et al 2006). 'SDSS' or 'UVES' indicates we 
have obtained w, 11196 from Sloan Digital Sky Survey spectra or from VLT/UVES spectra. 
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Figure 1. Correlation between absorber metallicity, [M/H], and Mg n rest-frame equivalent width, W; 



. The upper velocity scale is the approximate velocity 
spread calculated simply from Av [kms -1 ] x 70 [kms -1 A -1 ] X W; 12796 [A]. Left: Original data from Table[T]with additional [M/H]- and W? 2796 -uncertainties 
of 0.05 dex and 0.03 A added in quadrature (see text). The correlation is significant at the 4.2 tr level. Right: The [M/H] values have been corrected for 
cosmological metallicity evolution according to equation (TJ. This reduces the scatter, improving the correlation to the 4.7 <x level. 



Table 2. Statistics for different sub-samples. N- d is the sample size. [W; j and ([M/H]} are the median rest-frame Mgn equivalent width and the median 
metallicity, respectively. Columns 5-8 show the correlation statistics for the original data from Table[T](with increased errors) while columns 9-12 show the 
same statistics after the metallicities are corrected for redshift evolution (see text). P(r) is Kendall's t two-sided probability that no correlation exists between 
[M/H] and W^' 2796 . We convert this to a significance, 5(r), assuming Gaussian statistics. In all cases a Spearman test gives a more significant correlation. A 
and B are the intercept and slope, respectively, of the fit to the [M/H] versus Av data (see Section l3~TV 



Original data Redshift-corrected metallicities 



Sample 


AT. 


(wf 796 ) 


<[M/H]> 


P(T) 


S(t) 


A 


B 


P(t) 


S(t) 


A 


B 








[A] 


[dex] 


[l(r 3 ] 


[<r] 


[dex] 


[dex] 


[10- 3 ] 


[<r] 


[dex] 


[dex] 





Fiducial 


49 


1.38 


-0.99 


0.026 


4.2 


-4.29 ± 0.38 


1.69 + 0.20 


0.0032 


4.7 


-3.24 ± 0.35 


1.62 ±0.18 


1 


[Zn/H] only 


36 


1.33 


-0.97 


0.48 


3.5 


-4.28 ± 0.46 


1.70 + 0.24 


0.038 


4.1 


-3.14 ±0.40 


1.60 ±0.21 


2 


DLAs only 


40 


1.42 


-1.02 


0.028 


4.2 


-3.85 + 0.28 


1.41 ±0.14 


0.0011 


4.9 


-2.86 ± 0.28 


1.38 ±0.14 


3 


W; ,27% > 0.7 A 


42 


1.69 


-0.93 


2.3 


3.1 


-4.96 ± 0.50 


1.99 + 0.25 


1.7 


3.1 


-3.70 ± 0.44 


1.83 ±0.22 


4 


1 + 2 + 3 


26 


1.57 


-0.97 


2.2 


3.1 


-4.25 + 0.35 


1.59 + 0.17 


0.35 


3.6 


-2.97 ± 0.34 


1.44 ±0.16 


5 


W? 2196 < 1.4 A 


25 


0.96 


-1.20 


9.7 


2.6 


-4.54 ± 0.70 


1.91 ±0.39 


1.8 


3.1 


-3.67 ± 0.63 


1.94 ±0.35 


6 


W A2796 > ! 4 A 


24 


2.21 


-0.73 


100 


1.6 


-5.88 ± 1.32 


2.33 + 0.58 


17 


2.4 


-5.33 ±0.88 


2.50 ± 0.39 


7 


Zabs < 1-4 


24 


1.56 


-0.82 


20 


2.3 


-3.98 + 0.57 


1.64 ±0.30 


46 


2.0 


-3.28 ± 0.60 


1.67 ±0.31 


8 


Zabs > 1-4 


25 


1.33 


-1.19 


0.31 


3.6 


-4.16 + 0.36 


1.52 ±0.19 


0.048 


4.1 


-3.09 ± 0.34 


1.53 ±0.17 


9 


No H 2 


43 


1.29 


-1.00 


0.11 


3.9 


-4.52 ± 0.46 


1.83 ±0.24 


0.025 


4.2 


-3.38 ±0.42 


1.71 ±0.22 


10: NoH 2 or21-cm 


35 


1.25 


-0.99 


0.92 


3.3 


-4.80 ± 0.54 


1.97 ±0.29 


0.40 


3.5 


-3.49 ± 0.48 


1.78 ±0.25 


1 1 : H2 systems 


6 


2.78 


-0.84 


15 


2.4 


-4.06 ± 0.63 


1.44 ±0.29 


91 


1.7 


-3.26 ± 0.75 


1.55 ±0.32 


12: 21-cm systems 


9 


1.38 


-1.20 


95 


1.7 


-3.65 ± 0.68 


1.37 ±0.33 


12 


2.5 


-2.74 ± 0.63 


1.30 ±0.31 



sample, may somewhat alter its slope. To explore whether this is in- 
deed the case, we make a simple correction to the metallicities in 
Table[T]according to 

[M/H] corr = 0.52 + [M/H] + 0.27 z abs . (1) 

The numerical coefficients are the y-intercept and slope from an 
unweighted least squares fit to [M/H] versus z abs using the fidu- 
cial sample. This has the effect of correcting each metallicity by 
the expected mean value at the absorber's redshift. The slope and 
intercept are consistent with those derived from larger samples of 
DLAs covering simil ar redshift ranges (e.g. Proc haska et al.l 2003a ; 
Kulka rnTet al.l I2005T) . Fig. QJright) shows the redshift-corrected 
metallicities versus W^ 2796 and Table|2]gives the redshift-corrected 



statistics: the general scatter around the [M/H]-W/ 2796 correlation 
clearly decreases upon correction; the table shows that the signifi- 
cance of the [M/H]-W; 127% correlation increases in almost all sub- 
samples once the metallicities are corrected for redshift evolution. 

Visually, the correlation in Fig. [TJleft) appears tighter for the 
lower- W; 12796 half of the sample compared to that for the upper 
half. That is, the variance in [M/H] about the general [MfH]-W^ 2796 
trend appears smaller at low-W; ,27% . However, an F-test reveals 
only marginal evidence for this: a fit of [M/H] versus W; 12796 is per- 
formed [see equation {2); Section [3TT| , the sample is split into two 
sub-samples at W^ 2796 = 1.4 A and the variance in [M/H] about 
the fit is calculated for each sub-sample. The ratio of the vari- 
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ances about the fit is 1.83, but this should occur 15 per cent of the 
time by chance alone given the sizes of the sub-samples. Using 
the redshift-corrected metallicities [Fig. [Uright)], the ratio of the 
variances reduces to 1.64, which has an associated probability of 
24 per cent. The marginally increased scatter at high W; 12796 does 
reduce the significance of the correlation in the high-W; 12796 sub- 
sample compared to the low-lV; 12796 sub-sample (see Table|2j; one 
possibility is that Fig. [TJ particularly the redshift-corrected version 
in Fig. [TJ right), defines a lower bound 'envelope' to the metallicity 
at a given W^ 2196 instead of a normal correlation. This should of 
course be combined with an upper bound envelope at high metal- 
licities since arbitrarily high values are unphysical. A larger sample 
with robust control over metallicity uncertainties would be required 
to identify these features and to distinguish between a true correla- 
tion and a lower bound envelope. 

We have also split the sample into two redshift bins about the 
fiducial sample's median, z abs = 1-4. Table[2]shows that the [M/H]- 
W r ' 127% correlation is still well-defined in the high-z abs sub-sample 
but it is somewhat less statistically significant at low-z abs . However, 
the fits to both sub-samples, as described in Section [3~Tl yield con- 
sistent results. 

Three main objections might be raised to some of the data we 
use in Table[JJ (i) Some of the metallicities are not derived from Zn 
and so might be susceptible to dust-depletion effects, even though 
we have made a crude correction to any metallicities derived us- 
ing the heavily depleted Fe and Cr ions; (ii) Absorbers with low 
VK r ' 127% might contain a significant number of velocity components 
which are not completely saturated. While Wf 2796 still provides a 
measure of velocity spread, Av, in these systems, W^ 2796 will also 
depend somewhat on the optical depths of the unsaturated compo- 
nents. Therefore, as the Mgll velocity components weaken at lower 
[M/H], one might expect a correlation between [M/H] and W^ 2796 
which has little to do with the correlation found at higher W r 12796 . 
One might consider that only systems with W? 2796 > 0.7 A should 
be used since the profile is even more likely to be completely sat- 
urated; (iii) Finally, the reliability of sub-DLA metallicities is not 
as well established in the literature as that for DLAs. It is possible 
that we do not fully understand the ionization corrections required 
in some cases or, indeed, the magnitude of the systematics involved 
in deriving A'(Hi) in sub-DLAs. For example, Voigt profile 7V(Hi) 
estimates closer to the flat part of the curve of growth naturally have 
larger random uncertainties and may become increasingly, system- 
atically influenced by weaker Ly-ff forest lines. However, Table [2] 
shows that the [M/H]-lV r ' 12796 correlation is fairly robust against 
removal of non-Zn metallicities, systems with VV r ' 12796 < 0.7 A or 
sub-DLAs. Removing all three cases together also leaves a > 3cr 
correlation (> 3.5 cr after correcting the metallicities for redshift 
evolution). 

Note that our sample contains only one system with [M/H] < 
-1.7. This is partially a selection effect since the Znn lines become 
undetectably weak at low metallicities. However, very few DLAs 
and sub-DLAs are known with lower metallicities at redshifts 0.2 < 
Zabs < 2.6 and comparison with other more complete samples (e.g. 
iKulkarni et al.ll2005l:lAkerman et al.ll2005l) shows that our sample 
is fairly representative of the overall DLA metallicity distribution. 
Finally, we have included 3 systems in Table [JJ which are within 
5000 km s" 1 of the emission redshift of the background QSO. These 
systems do not show any signs of being associated with the QSO 
central engine and so are included here. Removing them from the 
sample makes a negligible difference to our results. 



3 DISCUSSION 

The [M/H]-W; 12796 correlation found in Fig.[JJis a direct demonstra- 
tion that the metallicity of strong Mg n absorbers is closely related 
to their kinematics (but not necessarily the kinematics of the host 
galaxy; see below). Since the Mgn A2796 transition is so strong, 
most velocity components in absorbers with W? 2796 > 0.3 A are 
saturated and so VV; 12796 measures the total velocity spread among 
the components, Av. Equivalently, it is a mea sure of the number of 
veloc ity components across the profile (e.g. IPetitiean & Bergeronl 
Il990h - Thus, the [MJH]-Wf 796 correlation demonstrates that the 
mechanism responsible for producing individual Mgn velocity 
components and for dispersing them over velocity ranges Av ~ 20- 
1000 km s~ is related to - and possibly also determines - the metal- 
licity of the absorber. 

The results in Fig.QJb ear on the apparent discrepancy between 
the conclusions drawn bv lNestor et alj d2003l). whose results sug- 



gest a [M/H]-IV/ 2796 correlation. and lYork etallfeOOot) whose anal- 
ysis suggests an anti-correlation. Our observed [M/H]-^ 12796 cor- 
relation supports the former and see ms inconsist ent with the latter. 
One possible explanation may be that lYork et al.1 assume a constant 
dust-to-gas ratio as a function of A'(Hi), i.e. E(B — V) oc A'(Hi). 
They base this assumption on their observation that the Mg n ab- 
sorbers they study seem to redden the background QSOs with, 
on average, an SMC-like dust extinction law. Even if this proves 
to be true for strong Mgn absorbers, this relation has certainly 
not been demonstrated in DLAs or sub-DLAs. Indeed, from the 
detection of SMC-like dust-reddening in > 300 SDSS DLAs at 
z a bs > 2.2, Murphy et al. (2007, in preparation) find that the 
mean E(B - V) ~ 7 x 10~ 3 appears not to change over the range 
20.3 < log[A'(Hi)/cm -2 ] < 2 1.7. One cav eat may be important 
here: possible magnitude bias. iNestor et alj assumed that 7V(Hi) 
is not a strong function of W; 12796 based on the observations of 
lRaoetalj |2006). The data in Table[JJalso reveal no apparent trend 
in /V(Hi) with W,' 12796 . However. lRao et all selected relatively bright 
background QSOs which could provide reasonable signal-to-noise 
rati os in ultrav iolet (UV) Hubble Space Telescope (HST) spec- 
tra. |%rkeUil]_suggest that dusty absorbers are therefore unlikely 



to exist in 



Rao et alj 's sample and that th is dust-bias (relative to 



deeper surveys like SDSS) may invalidate Nes tor et alj 's assump- 
tion. From this point of view, it may not be sur prising that our re- 
sults agree with the conclusion of INestor et alj since most, but not 
all, of the absorbers in our sample occult relatively bright QSOs. 



3.1 Comparison with higher redshift DLAs 

At higher redshifts (1.6 < z abs < 3.0). IWolfe & Prochaskj jl998l) 
tentatively noted a co rrelation between A v and metallicity from a 
sample of 17 DLAs. IPeroux et al] d2003h found a similarly ten- 
tative [Zn/H]-Av correlation from a larger sample of DLAs and 
sub-DLAs over the redshift range 1.4 < z abs < 4.5. However, 
lLedouxetaH l l2006l hereafter lL061) very recently found a relatively 
tight and clear [M/H]-Av correlation from a sample of UVES spec- 
tra containing 70 DLAs and sub-DLAs over the range 1.7 < z abs < 
4.3. They defined Av from different transitions of several low ion- 
ization species (Oi, Sin, Fen, Crn, Sn) with moderate optical 
depths such that the strongest velocity component absorbed be- 
tween 10 an d 60 per cent of the continuum. Since many of the 
DLAs in the IL06I sample are at z abs > 2.6 where the Mg n lines 
are not observable with optical spectroscopy, a complete compari- 
son between the sample and ours is difficult. Nevertheless, it 
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is clear that, since W; i27% is a measure of Av, the [M/H]-W; 12796 
correlation in Fig.Q]persists at higher redshifts. 

|L06| fit their [M/H]-Av data using the least squares bisector 
method of llsobe etai] (fl990). finding the best-fit relationship to be 
[M/H] = (1.55 ± 0.12)logAv[kms-'] - (4.33 ± 0.23). For com- 
parison, we performed fits using the same bisector technique after 
converting our Wf 2796 values to Av according to the simple relation 
Av[kms-'] » 70[kms _1 A" 1 ] x Wf 96 [k]. The proportionality 
constant is derived from a fit of W; 12796 versus Av, with the inverse 
of the variances on W r ' 127% used as weights, for the 18 absorbers 
which appear in both our sample and that of lLOd The fiducial sam- 
ple is best-fitted by the relation 



[M/H] = (1.69 + 0.20) log Av [kms~'] - (4.29 ± 0.38) , 



(2) 



which is consistent with the relation found bv lLOq . 

Table [2] shows the intercept A and slope B of the [M/H]-Av 
relationship for the different sub-samples. In particular, we find no 
evidence for significant differences in A or B in t he lo w- and high- 
Zabs sub-samples. This is similar to the results of L06 who find the 
slope and intercept to be statistically consistent in two redshift bins 
split at z abs = 2.43. We also note that the median W? 2796 and [M/H] 
of our low-z a b s sub-sample are marginally higher than at high-z a (, s : 
(Wf 796 ) = 1.56 A and <[M/H]> = -0.82 for the former compared 

with (W? 2796 ) = 1.33Aand<[M/H]> = -1.19 for the latter. Again, 
L06 find similar trends. 



3.2 Mass-metallicity correlation or anti-correlation? 

The [M/H]-W? 2196 (or [M/H]-Av) correlation in DLAs and sub- 
DLAs provides an important link between the kinematics of the 
absorber and the metal-enrichment history of the absorber. How- 
ever, one would like to go a step further and link the metallicity 
with the mass of the halo in which the absorber resides. This step is 
potentially confusing. For example, many authors have, in the past, 
tacitly assumed that the absorption-line kinematics are a reliable 
tracer of the host-galaxy kinematics. That is, Av is assumed to pos- 
itively correlate with the circular velocity of the host-galaxy. By 
presuming this detail of the absorber-galaxy connection - which 
should instead be observationally determined - one would conclude 
that W; 127% (and Av) is correlated with the galaxy mass, make 
this latter assumption and use it to constrain the implied mass- 
metallicity relationship. Indeed, they find this to be consistent with 
the luminosity-metallicity relationship derived from local galaxies. 

However, the assumption that galaxy or halo mass is corre- 
lated w ith W? 279 6 is c hallenged by new observations. iBouche et al.l 
d2006l hereafter lBOil identified 1806 Mgn systems with Wf 796 > 
0.3 A and ~ 250000 luminous red galaxies (LRGs) within projected 
co-moving distances of 13/T 1 Mpc of the absorbers. The ratio of the 
absorber-LRG cross-correlation to the LRG-LRG auto-correlation 
provides a measure of the bias-ratio between absorbers and LRGs. 
In hierarchical structure formation scenarios the bias scales with 
halo mass, so the clustering of LRGs around the Mg n absorbers 
prov ides a measure of the absorber halo mass. With this approach, 
[Bod found that halo mass and W r >2796 are anti-correlated. They also 
found that by combining the previously observed luminosity- and 
W'* 2796 -dependence of the absorber cross-section, or by combin- 
ing the cross-section's luminosity-dependence with the observed 
incidence probability of strong Mgn absorbers, d 2 N/dz il bs < i^ t2796 , 
one finds that luminosity (or m ass) and W r ' 12796 should be anti- 
correlated. IProchter et alj |2006) reached a similar conclusion by 



interpreting the observed redshift evolution of the number density 
of strong Mg n systems with simple cross-secti on arg uments. 

Of course, the main difference between the|B06J absorber sam- 
ple and the one studied here is that B06 used strong Mg n absorbers 
while we study DLAs and strong sub-DLAs. It is possible that a 
mass-W^ 12796 correlation does exist for DLAs if the sharp drop in 
halo-mass observed bv lB06l for Mgn systems with VV r ' 127% > 1.0 A 
were dominated by systems with neutral hydrogen column densi- 
ties below the lower limit imposed here, i.e. N(H i) = 3 x 10 19 cnr 2 . 
However, this seems very unlikely given that the fraction of strong 
Mgn systems which are DLAs [i.e. have /V(Hi) = 2 x 10 20 cnr 2 ] 
increa ses from ~ 15 to ~ 65 per cent over the range Wf 2796 = 0.6- 
3.3 A (|Rao et alj20061) . It is therefore important to confirm or refute 
the |B06j results, which derive from SDSS Data Release 3. Bouche 
et al. (2007, in preparation) will present an analysis based on a sig- 
nificantly increased sample of absorbers from Data Release 5. 



3.3 Are strong Mg n absorbers caused by outflows? 



[Bod interpret their observed mass-W; 12796 anti-correlation as evi- 
dence that outflows become the dominant mechanism for produc- 
ing Mgn absorbers as W? 2796 increases beyond ~ 1.0 A. They rea- 
son that outflows are more easily ejected from lower mass star- 
forming galaxies and so, when viewed along a random QSO sight- 
line, are more likely to produce more Mgn velocity components 
over a larger velocity range. In this outflow picture, W; 12796 and 
the number of Mg n velocity components may be inversely related 
to the mass and could depend mor e stro ngly on the recent star- 
formation history of the host-galaxy, mention that some frac- 
tion of DLAs could be caused by outflows, but they again attribute 
the outflow velocity spread to be directly related to the depth of the 
host galaxy's potential well. 

An initial objection to the outflow model for DLAs and strong 
Mg ii systems might be that such hot, highly ionized material may 
not produce velocity components of cold gas which can be traced 
with low-ionization species such as Mgn. However, the outflows 
from local starburst g alaxies are well-known to c ontain clouds 
of cold gas and dus t dLehnert & Heckman[ 19961 ; Martini 1 19991 : 
lHeckman etai]|200d : iRupke et all 12005k iMartinf l2006h . Dense, 
molecular gas is also known t o be entrained in smaller scale out- 
flows (e.g. Nak ai et ail 1 19871) and even the ionization cones of 
nearby active galactic nuclei which s how starburs t activity (e.g. 
Irwin & Sofue 1992; Curran et al. 1999). Moreover, iNorman et al.l 
( 1996) identified Mg n absorption in spectra of two QSO sight-lines 
with impact parameters 25 h~ [ kpc and 55 h~ [ kpc near the galaxy 
NGC 520, a local starburst with super-winds, a disturbed morphol- 
ogy and filamentary Ha emission along its minor axis. Thus, it 
seems clear that outflows from star-forming galaxies do contain 
cold gas. 

If indeed many high-W; 12796 Mg n absorbers are caused by out- 
flows and these outflows do indeed entrain cold clouds (or cold 
gas clumps cool out of the hotter outflow medium), do we ob- 
serve tracers of the cold, dusty gas, such as molecular hydrogen? 
H 2 is identified via UV absorption in the Lyman and Werner bands 
and is known to ar ise in very cold, dusty velocity components (e.g. 
Led oux et ai] |2003h . It is therefore notable that 3 of the 5 highest- 
W^ 12796 systems in Fig.[TJare known to contain H 2 . Indeed, the two 
highest- W- 12796 systems those towards Q 0013-0029 and Q0551- 
366, both contain H 2 dLedoux et ai] |2002t IPetitiean et al] [2002) 
and have Wf- 196 values at the extreme of the Wf 196 distribution 
obtained from large SDSS Mgn surveys (e.g. iNestor et alj|2003t 
IProchter et ai1l2006h . The median W? 2796 for the H 2 absorbers in 
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Figure 2. Montage of the velocity structures of strong transitions in the 13 known tb -bearing DLAs/sub-DLAs in the literature (6 are in our sample) and the 
21-cm absorbers in our sample. The transition used to illustrate the extent and distribution of absorption components is marked on each panel. In most cases 
the Fe n /12382 line is used. Fen ,42600 is used for the absorber towards Q 0458-020. When Fen /12382 was not available (particularly at high-z a (, s ), Oi 302, 
Sin A 1526, Alii /tl670 or Mgn /12796 is used. All systems are shown on a common velocity scale registered to the redshift of the strongest H2 or 21-cm 
component in each system. The black histogram is the data normalized by a local continuum fit while the grey/red line near zero flux shows the normalized 
1 o" error array. The approximate extent of the transitions is marked by the bracket below the zero level. Note the similarity in the profiles for many of the 
H2-bearing systems. These broad and complex structures may be characteristic of outflows. The 21-cm profiles generally do not show such outflow signatures, 
but the sample is too small for strong conclusions to be made. 



our sample is 2.8 A compared with 1 .4 A for the sample as a whole 
(Tabled. Even with only six H 2 systems, a Kolmogorov-Smirnov 
test gives a probability of just 9 per cent that the equivalent widths 
of the H2 absorbers are drawn from the same parent population as 
the rest of the sample. 

In Fig. [2] we compare the velocity structures of all the 13 
known H2 -bearing DLAs and sub-DLAs. It is striking how similar 
the velocity structures are for the systems towards Q 05 15-4414, 
Q 055 1-366 and Q 0013-0029: each is spread over Av > 
700kms~' and comprises at least three groups of stronger veloc- 
ity components separated by few (if any) weaker components. The 
systems towards Q 1444+014 and Q2343+125 have very simi- 
lar profiles and these are not dissimilar to the above three. Also, 
the systems towards Q 2348-011 and Q 0528-250 both have two 
broad, saturated regions separated by a narrow region of low opti- 
cal depth. Such st ruct ures and symmetries have been discussed by 
Bond et alj J200lh and lEllison etaU b003h as being the character- 
istic signs of absorption in galactic outflows. It is also conceivable 
that th ese profiles are pr oduced in galaxy mergers (but see discus- 
sion in lBond et ai]|200ll) . 

Although no strong conclusions can be reached due to the 
sample size, Fig.[2]nevertheless suggests that approximately half of 
the known H 2 -bearing systems may arise in outflows. If indeed cold 
gas and dust is entrained in outflows and the gas is well-shielded 
then H2 should form on the surface of dust-grains. Even if this hap- 
pens in very few clouds, the large number of clouds across these 
profiles significantly increases the probability of observing one 
or more H2-bearing clouds towards any given background QSO. 



Hence, even if these outflow-like absorption profiles are themselves 
rare, H2 may be found within them in a large fraction of cases. 

Fig. [2] also shows the velocity structures of the 21-cm ab- 
sorbers in our sample for which UVES spectra were available. 
These velocity structures are typical of those in the general popula- 
tion, seemingly consistent with the fact that the 21-cm absorbers are 
generally indistinct from most other systems in the [M/fTJ-W; 12796 
plane (Fig. [T). Like H 2 absorption, 21-cm absorption should arise 
in cold clouds, but there is no evidence in Fig.|2]to suggest that 21- 
cm absorbers also arise in the complicated and broad outflow-like 
profiles like the H 2 systems. However, the 21-cm sample is very 
small and is dominated by low-z a b s systems; outflows are likely to 
be less common below z a b s ~ 1 due to the general decline in the 
global star-formation rate density (see, e.g jBond et aL 2001). 

IL06| note that, under the assumption of a mass-metallicity re- 
lationship, the absorbers with the highest equivalent widths and 
metallicities will be most easily detected in direct imaging obser- 
vations. This prediction is also natural in the outflow picture where 
Wf 1196 and the number of Mgn velocity components are related 
to both the mass and the recent star-formation history of the ab- 
sorber host-galaxy. While one expects the highest Wf 2796 and high- 
est [M/H] systems to be hosted by lower mass galaxies, they should 
be actively star-forming and show strong Ha emission. If indeed 
H 2 -bearing systems more often probe star-forming galaxies with 
outflows, they may be more reliable tracers of the cosmological 
evolution of metallicity than other DLAs, many of which might 
not be closely linked to star-formation sites/process es in high-z 
galaxi es. We have already advan ced this possibility in lCurran et al.l 
j2004l) and lMurphv et al j (2004) based on a comparison of H2 ab- 
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sorber metallicities with those of the general DLA population. The 
H 2 absorbers showed a faster, more well-defined increase in [M/H] 
with decreasing redshift, although the sample was quite small. A 
much larger sample of H 2 -bearing systems is required to compute 
the 7V(H i)-weighted metallicity evolution in a way comparable with 
general DLA samples. 



4 CONCLUSIONS 

For 49 DLAs and sub-DLAs with published metallicities we have 
measured or gathered from the literature their Mgn rest-frame 
equivalent widths to search for an empirical relationship between 
the kinematic spread, Av, in the absorbing metal-line velocity com- 
ponents and the total absorber metallicity. The vast majority of 
DLAs and strong sub-DLAs have Wf 1196 > 0.3 A which means 
that most Mg n velocity components across the absorption profile 
are saturated. Thus, W? 2796 provides a simple measure of Av which 
can be derived even from relatively low-resolution QSO spectra. 

We find a correlation between [M/H] and W? 2796 at the signifi- 
cance level of 4.2 cr [Fig.[TJleft)] which increases to 4.7 a when the 
[M/H] values are corrected for mild cosmological evolution over 
the sample's redshift range of 0.2 < z a u < 2.6 [Fig.QJright)]. Even 
for our most conservative sub-sample, where we include only 26 
DLAs (no sub-DLAs) with Zn-metallicities and W? 2796 > 0.7 A, 
we still find a correlation at > 3 a (> 3.5 <x after redshift-correction 
of [M/H]). After converting the W; 127% values to Av, the slope of 
our [M/H ]-log Av relat i onship is consistent with that reported re- 
cently by Ledo uxet al.l J2006h for generally higher-z a b s DLAs and 
sub-DLAs, suggesting that the [M/tLJ-lVJ 12796 correlation persists 
up to z a bs = 4.7. Although we find a lower statistical significance 
for the correlation below z abs = 1 .4, the fitted slope and intercept of 
the [M/H]-log Av data is similar for the z abs < 1 .4 and z a bs > 1 -4 
sub-samples. Finally, we find marginal (~ 1.4 a) evidence for in- 
creased variance in metallicity above WA 12796 = 1.4 A. However, the 
evidence weakens once the [M/H] values are corrected for evolu- 
tion; other non-intrinsic sources of metallicity scatter may need to 
be addressed before such an effect can be verified. 

If one assumes that the metal-line absorption kinematics trace 
the host-galaxy kinematics - that Av is proportional to the host- 
galaxy circular velocity - then it is natural to translate the ob- 
served [M/HJ-W/ 2796 cor relation into a corre lation between [M/H] 
and galaxy mass. Indeed, iLedoux et alj fc006t) make this assump- 
tion and find that the resulting [M/H]-mass correlation is consistent 
with the [M/H]-luminosity relationship derived from local galax- 
ies. However, by studying th e clustering of galaxies around strong 
Mg ii absorbers in the SDSS. lBouche et alj (2006) find that mass is 
anti-correlated with IV/ 2796 . This result, together with the [M/H]- 
W r ' 12796 relationship observed here, implies that the absorber metal- 
licity is anti-correlated with halo or galaxy mass. It is th erefore 
important to confirm or r efute the resu l ts of Bouche et al.1 and to 
address the assumption of Ledo ux et al.1 d2006 ) with observations. 

Finally, we note that about half of the known H 2 absorbers 
have very broad velocity structures which show distinct distribu- 
tions of velocity components (Fig. [2). The patterns, groups and 
symmetries evident in these cases are consistent with an outflow 
origin for the bulk of the components. Since outflows can entrain 
cold, dusty gas and since they should provide a large number of 
velocity components, the probability of finding one or more H 2 - 
bearing components in any given outflow-driven absorber might 
be quite high, even though H 2 -bearing components are themselves 
quite rare. This suggests the targeting of high-W; 12796 absorbers to 



identify more H 2 -bearing systems. Given the [M/HJ-W^ 2796 corre- 
lation we o bserve, this may be sim ilar to the metallicity-selection 
approach of Petitjean et al 1 d2006h . However, it has the advantage 
that one can determine W;' 2796 with only moderate-resolution spec- 
tra - or infer it from the equivalent width of bluer transitions when 
z a b s <; 2.5 - which need not be blue- or UV-sensitive (e.g. SDSS 
spectra) . One can then follow-up the highest IV; 12796 systems with 
high-resolution, UV- or blue- sensitive spectra to search for H 2 ab- 
sorption. 
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